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ABSTRACT 

We present an in-depth study of the High Mass X-ray Binary (HMXB) GX 301-2 during its pre-periastron flare using data from the 
XMM-Newton satellite. The energy spectrum shows a power law continuum absorbed by a large equivalent hydrogen column on the 
order of 10 24 CUT 2 and a prominent Fe Ka fluorescent emission line. Besides the Fe K a line, evidence for Fe K/J, Ni Ka, Ni KJ3, 
S Ka. Ar Ka, Ca Ka, and Cr Ka fluorescent lines is found. The observed line strengths are consistent with fluorescence in a cold 
absorber. This is the first time that Cr Ka is seen in emission in the X-ray spectrum of a HMXB. In addition to the modulation by 
the strong pulse period of -685 sec the source is highly variable and shows different states of activity. We perform time-resolved as 
well as pulse-to-pulse resolved spectroscopy to investigate differences between these states of activity. We find that fluorescent line 
fluxes are strongly variable and generally follow the overall flux. The N H value is variable by a factor of 2, but not correlated to 
continuum normalization. We find an interval of low flux in the light curve in which the pulsations cease almost completely, without 
any indication of an increasing absorption column. We investigate this dip in detail and argue that it is most likely that during the 
dip the accretion ceased and the afterglow of the fluorescent iron accounted for the main portion of the X-ray flux. A similar dip was 
found earlier in RXTE data, and we compare our findings to these results. 

Key words, stars: neutron (GX 301-2) - X-rays: binaries - Accretion, accretion disks 


1. Introduction 

An X-ray flare at the coordinates of the High Mass X-ray Binary 
(HMXB) GX 301-2 was first detected in 1969 during a balloon 
experiment (Lewin et al., 1971). During further balloon flights 
McClintock et al. (1971) found that the source was persistent. 
Since then GX 301—2 has been observed by most X-ray instru- 
ments and has stirred a lot of interest in the community as it 
shows a strongly variable flux and a very large absorption col- 
umn. The source shows an erratically changing pulse period of 
~ 685 sec, detected by White et al. (1976) and in the last few 
years monitored by the Ferm//GBM Pulsar Project (Finger et al., 
2009) 1 , clearly identifying the compact object as a neutron star. 
The system has an orbital period of P or b = 41 ,506±0.003 d at the 
reference time 7’o(MJD) = 43906. 06±0.1 1 and the orbital period 
is changing with P (irb = (-3.7 ± 0.5) x 1 0 -6 s s" 1 (Doroshenko 
et al., 2010). These values were determined by accurate pulse 
timing and are in agreement with earlier results from Koh et al. 
(1997) who used the regular X-ray outburst of GX 301-2. The 
eccentricity of the orbit was determined to be e = 0.47. 

The neutron star accretes matter from the optical companion, 
which was identified as the B-type giant Wray 977 (Vidal, 1973). 


1 see also http://www.batse.msfc.nasa.gov/gbm/science/ 

pulsars/lightcurves/gx361m2 . html 


The exact type of the companion and its distance was long dis- 
puted. but the most recent measurements by Kaper et al. (2006) 
indicate a distance of 3 kpc and identify Wray 977 as a B I Ia+ hy- 
pergiant with a stellar-radius of 62 R Q and a mass of 43 ± 10M e . 
We adopt these values throughout the paper, as X-ray accretion 
models by Leahy & Kostka (2008) also seem to indicate similar 
values. 

The regular and bright X-ray flares of GX 301-2 were found 
to occur - 1 .4d before the periastron passage, i.e., before the neu- 
tron star reaches the densest part of the stellar wind. Standard, 
symmetric wind accretion models fail to explain pre-periastron 
flare. Pravdo & Ghosh (2001 ) suggested that a circumstellar disk 
is formed around Wray 977, from which the neutron star ac- 
cretes additional matter when it passes through. If the disk is not 
aligned with the orbital plane of the neutron star, a pre-periastron 
flare could be observed. More recent calculations (Leahy , 2002; 
Leahy & Kostka, 2008) show that it is also possible that the neu- 
tron star is trailed by an accretion stream, which it overtakes 
shortly before periastron. This model self-consistently explains 
the observed dependence of flux and column density on the or- 
bital phase. It is the regular pre-periastron flare which distin- 
guishes GX 301-2 from most other HMXB. During the flare 
the X-ray flux increases by up to a factor of 25 (Pravdo et al., 
1995; Rothschild & Soong, 1987), reaching luminosities of up 
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Fig. 1. Sketch of the orbit of GX 301-2, using the parameters given 
by Koh et al. (1997). The green line shows the average RXTE - ASM 
light curve as function of orbital phase. The phase of the XMM-Newton 
observation discussed here is highlighted in red. 


to 1 Crab in the 15-50keV band and some lOOmCrab in the 
softer 2-10keV band. X-rays with energies below ~5keV are 
less variable over the orbit as they are strongly affected by scat- 
tering and absorption due to the high column density (Leahy & 
Kostka,2008). 

Figure 1 shows the orbit of GX 301-2 and the RXTE - ASM 
light curve folded onto an orbital period of 41 .5 d using the 
ephemeris by Doroshenko et al. (2010). Clearly visible is the 
pre-periastron flare, peaking around phase tp = 0.91 . The plot 
assumes a size of 62 R 0 for the optical companion (Kaper et al., 
2006) . 

The broad-band X-ray spectrum of GX 301-2 can be de- 
scribed by a power law with a high-energy cutoff around 20 keV 
and a cyclotron resonant scattering feature (CRSF) around 
37 keV (Endo et al., 2002). Kreykenbohm et al. (2004) showed 
that the spectral parameters, e.g., the energy of the CRSF, vary 
with pulse phase, and therefore concluded that the X-rays are 
most likely produced in an accretion mound close to the neutron 
star surface. In the softer X-rays a very strong Fe K a fluorescent 
emission line is visible together with a photo-electric absorp- 
tion with an equivalent hydrogen column 1% of up to 10 24 cm"' 2 
(Leahy et al., 1988, 1989). Using Chandra data, Watanabe et al. 
(2003) detected a fully resolved Compton shoulder at the low- 
energy flank of the Fe Ktr line. These authors could model the 
shoulder by simulating the Compton scattering of photons in a 
dense, cold gas-stream. 

The slow pulse period of GX 301-2 is strongly variable 
over long timescales (Pravdo & Ghosh, 2001; Evangelista et al., 
2010). It was first detected with P ~ 700sec, but changed in 
the last 40 years to values as low as 675 sec and is today around 
P ~ 685 sec. Although the period is variable, GX 301-2 usu- 
ally show's clear pulsations throughout its orbit. Recently, how- 
ever, GogU§ et al. (201 1), discovered a dip in the RXTE X-ray- 
light curve during which GX 301-2 ceased to show pulsations. 
The data showed that exactly one pulse was missing and that the 
pulsed fraction was gradually reduced before and gradually in- 
creased after the dip. Interestingly, the spectrum softened signif- 


icantly during that dip, contradicting the common picture that a 
very dense blob, absorbing all soft X-rays, was passing through 
the line of sight. Gogii§ et al. (201 1) speculate that due to the 
cessation of accretion the black body radiation from the polar 
cap became visible, effectively softening the spectrum and re- 
ducing the pulsed fraction, as the polar cap is larger and un- 
dergoes stronger gravitational light bending than the accretion 
column. Similar dips have already been observed in the light 
curves of Vela X- 1 (Kreykenbohm et al ., 2008) , but never before 
in GX 301-2. 

The strong absorption during the bright pre-periastron flare 
makes GX 301-2 an ideal object to study the X-ray producing 
and absorbing regions around a neutron star, especially focusing 
on the many fluorescent lines. XMM-Newton is the ideal instru- 
ment for this observation, with its large effective area and good 
spectral resolution. We obtained a ~47ksec XMM-Newton ob- 
servation of which the results are presented here. 

In the next Section we will describe the data used. Section 3 
will cover the timing analysis and in Sect. 4 the different spectral 
analyses will be discussed. In Sect. 5 the results are discussed 
and a conclusion and an outlook are given. 

2. Data & Observations 

We present data taken with the XMM-Newton observa- 
tory on July 12th, 2009 (MJD 55024.103-55024.643, ObsID 
0555200401) at orbital phase ~0.91 of GX 301-2, i.e., during 
the pre-periastron flare. The exposure was 46.6 ksec, covering 
only a small part of the flare (see Fig. 1). The X-ray flux was 
strongly variable during the observation, with values between 
3.67 x 10~ 10 erg sec -1 cm" 2 and 2.47 x 10~ 9 erg sec" 1 cm" 2 in the 
2-1 0 keV energy band. Due to the expected high flux, the EPIC- 
MOS detectors were turned off to provide more telemetry for 
the EPIC-pn camera, which was operated in timing mode. In 
this mode the CCD is continuously read out, so that all spatial 
information in the read out direction is lost. In the other direc- 
tion only the central 64 columns are used. The timing mode al- 
lows the measurement of countrates up to 400 cts sec" 1 without 
detectable pile-up (Striider et al., 2001). As GX 301-2 shows 
very strong absorption no usable information can be gained from 
the Reflection Grating Spectrometers (RGS), which are mainly 
sensitive below 2.5 keV. We used the standard Science Analysis 
System (SAS) 10.0.0 to extract light curves and spectra. We 
tested for pile-up using the tool epatplot, but found no severe 
pile-up. The source data were extracted from EPIC-pn columns 
35 — 40 and background data were obtained from columns 10-17 
and 53-60. As high photon fluxes change the charge transfer in- 
efficiency (CTI) of the CCD detector, we also extracted spectra 
with the tool epfast, which takes these changes into account. 
We found, however, no significant differences to the standard 
pipeline. The data were barycentric and binary corrected, using 
the ephemeris provided by Doroshenko et al. (2010). 

3. Timing Analysis 

3.1. Light curve 

The top panel of Fig. 2 shows the 20 sec resolution light curve of 
GX 30 1 -2. The light curve shows clear flaring behavior, with a 
smaller flare in the beginning and a large flare starting ~30 ksec 
after the start of the observation. Clearly visible are also the sub- 
peaks induced by the pulse period of ~685sec. We performed 
a search for the pulse period via epoch folding and found a pe- 
riod of P = 685.0 ± 0. 1 sec, consistent with the monitoring data 


2 


F. Fiirst et al.: Fluorescent lines and absorption in GX 301-2 



i i i UK LLC 

10000 20000 30000 40000 

Time [sec] 


Fig.2. a) Lightcurve of GX 301-2 with 20 sec resolution in the 0.5- 
1 0.0 keV energy band. Vertical lines indicate the time regimes I— IV used 
in the later analysis, b) Hardness ratio between the continuum energy 
band C (8.5-1 0.0 keV) and the absorbed energy band A (4-.0-5.0keV). 
c) Hardness ratio between C and the Fe K a line band F (63-6.5 keV). 


from the Fermi/GBM Pulsar Project. In search for changes of the 
pulse period during the observation, we performed epoch fold- 
ing and phase connecting on parts of the light curve as short as 
1 5 ksec. We did not find a significant change of the pulse period, 
consistent with Ferm;/GBM data which shows that the pulse pe- 
riod was only slowly changing during that time. 

A close-up view of the light curve reveals that albeit the pul- 
sations are clearly visible, strong pulse-to-pulse variations are 
also present. The individual pulse shapes range from simple two- 
peaked profiles to more complex shapes, with multiple minor 
peaks superimposed on the general shape. An example of two 
different profiles is shown in Fig. 3 using light curves with 1 sec 
and 20 sec resolution. It becomes clear that the average pulse 
profile, also shown in Fig. 3, is highly complex. Additionally 
brightness changes independent of the pulses are visible in the 
light curve. These changes can be on the same order with respect 
to luminosity and duration as the variability due to the pulsation. 

Stronger and sudden changes in luminosity are also seen in 
the light curve and different luminosity levels are marked as re- 
gions I— IV in Fig. 2. These changes are also accompanied by 
changes in the pulse profile and the pulsed fraction. To investi- 
gate this in more detail, we folded each regime I— IV separately 
on the pulse period of P = 685.0sec and studied the normalized 
distribution of the countrate in each phase bin. As the overall 
countrate is different for every region we also normalized each 
light curve to its maximum countrate. The mean of the distribu- 
tion in each phase bin gives the value of the pulse profile in that 
bin. Figure 4 shows the results for all four regimes in a color- 
coded map. Region III shows the clearest pulse profile, where 
the mean countrate as well as the whole distribution is changing 
with phase. This behavior is different in region I. Here the pulse 
profile is also visible in the mean values but the values measured 
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Fig.3.Two examples of the variability of the pulse shape, left an exam- 
ple early in the observation, right one close to the end of the observation. 
Shown in black is the 1 sec light curve, in red a 20 sec light curve. The 
blue profile is the average pulse profile over the complete observation. 
Note that the y-axes are not the same scale. 



Phase Phase 

Fig. 4. Color-coded maps of the pulse phase distribution of normalized 
countrates in regions I-IV as indicated in Fig. 2. See text for details. 
Color represents the number of measurements for each countrate in the 
respective phase bin. Each phase bin was normalized to itself, so that 
they have all comparable color-codes. Colors are ranging from violet, 
corresponding to zero, to red, corresponding to one. The superimposed 
black line shows the mean value, i.e., the pulse profile. 


most often are only weakly phase dependent. The two peaks of 
the profile result from strong outliers and flares at the respective 
phases. In region II a similar effect is seen, but the overall pulsa- 
tions are even weaker. Region IV has a very short duration and 
therefore provides barely enough statistics for an interpretation. 
The behavior seems similar to the one seen in region III, but the 
pulse profile is not that well pronounced. 

Region II seems to be the region with the weakest pulsations 
and shows a very different behavior than region III. To investi- 
gate these differences further, we took a closer look at the 1 sec 
light curve during region II, as shown in Fig. 5. In the beginning 
the pulse is visible but it vanishes almost completely during the 
decline of the luminosity for ~3 rotations. Around t = 21 000 sec 
a sudden increase of the luminosity is accompanied by weak 
pulses but these disappear again after ~4 rotations of the neu- 
tron star, also during a dip in the brightness. Looking at the hard- 
ness ratio during that time (Fig. 2, lower two panels) allows us 
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Fig.5. Small part of the light curve between 0.3-12.0 keV, showing a 
vanishing of the pulsations. In black the light curve with 5 sec resolution 
is shown, in red with 40 sec resolution. Superimposed in blue is the 
average pulse profile for the whole observation. 


to deduce the spectral shape. The hardness ratio was calculated 
as HR(W, S) = (H-S)/(H + S), with //indicating the count rate 
in the harder and S in the softer energy band. When comparing 
the count rate in the 8.5-10keV continuum band C to the count 
rate in the strongly absorbed energy band A between 4—5 keV, 
no significant variability of the hardness HR(C,A) can be mea- 
sured during region II. Neither is a change compared to region 
I visible, indicating that the reduced flux is not due to increased 
absorption, which would primarily absorb the soft X-rays and 
would thereby make the spectrum overall harder. Contrary to 
that we see that the hardness HR(C, F) between the count rate 
in the continuum band and in the iron line energy band F be- 
tween 63-6.5 keV changes dramatically in region II. The strong 
drop of the hardness indicates that the flux of the iron line is not 
reduced as strongly as the continuum flux. 

3.2. Pulse profiles and pulsed fraction 

In the energy range of XMM-Newton , the average pulse profile 
is only weakly energy dependent (Fig. 6). At higher energies the 
smaller, secondary peak becomes relatively stronger compared 
to the first one, but still stays clearly weaker. Only at energies 
above 20keV the secondary peak gets almost as bright as the 
primary one (Kreykenbohm et al., 2004; La Barbera et al., 2005). 
This behavior is also confirmed in recent Suzaku data (Suchy 
et al., 201 1). 

Endo et al. (2002) and Tashiro et al. (1991) have shown that 
the pulsed fraction in the iron line energy band is significantly 
smaller than in other neighboring energy bands. Analysis of 
XMM-Newton data confirms this behavior, indicating that the 
matter producing the fluorescent iron line is distributed homo- 
geneously around the X-ray producing region. The pulsed frac- 
tion / was calculated using the formula 

_ max(P) - min(P) 
max(P) + min(P) 

with P denoting the pulse profile. The pulsed fraction is plotted 
in Fig. 7 for the four different regions defined above. Common 
to all four regions is the fact that the pulsations in the iron line 
band are notably reduced. In region IV the dip is least prominent, 
probably due to the overall low pulsations and the short duration 
of this region. 

Figure 7 also clearly shows that the pulsed fraction is vari- 
able with time, dropping significantly during region II, as already 
seen in Fig. 4. Some residual pulsed fraction is still measurable 



Phase 

Fig.6. Pulse profiles in different energy bands: a) 0.3-4.0keV, b) 4.0- 
6.2 keV, c) 6.2-6.6keV, and d) 6.6-15.0 keV. Each profile is shown 
twice for clarity. 
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Fig. 7. Pulsed fraction for different energy bands and for the 4 different 
regions I— IV. The dashed vertical line marks 6.4 keV, the energy of the 
Fe K a line, the dotted one 7.1 keV, the energy of the Fe K fi line. 


outside of the iron line band, but as the iron line band domi- 
nates the whole spectrum during that interval, almost no pulsa- 
tions are visible in the broad band light curve. The lowest energy 
bin, between 0.5^I.0keV is dominated by the strong absorption, 
smearing out the pulsations and reducing the signal-to-noise ra- 
tio (S/N). With increasing energy the pulsed fraction increases 
above 6.4 keV, with exception of the energy bin covering the 
Fe K/3 line at 7.1 keV. The Fe K/3 line is strong enough to con- 
tribute significantly to the flux in that bin, and the line is, as the 
Fe Ka line, only very weakly pulsed. 
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4. Spectral analysis 

4.1. Phase averaged spectrum 

When looking at the averaged X-ray spectrum of GX 301-2 
(Fig. 8a) one sees a large iron K a fluorescent emission line and 
an equally large absorption column, which reduces the countrate 
to effectively zero below 2 keV. The spectrum can be reasonably 
well modeled with a partially covered power law, with additional 
Gaussian lines at the Fe K a, Fe K /3, Ni Ka, Ca Ka, Ar Ka, and 
S Ka fluorescent line energies. We used a model of the form: 

(V® al x((dV H ,i+(l -c)N H o) x power law + Gaussian lines j (2) 

Where Nh,i and Nu ,2 are the two different parts of the partial 
coverer and 0 < c < 1 is the covering fraction. We applied the 
tbnew absorption model (Wilms et al., 201 1 , in prep.) 2 for both 
elements of the partial coverer. For the large columns measured 
in GX 301-2 the plasma becomes optically thick and some flux 
is scattered out of the line of sight. This effect is taken into ac- 
count with the cabs model, described by 

cabs(£) = exp(-/V H <x T (£)) (3) 

where o"t(£) is the Thomson cross section and N h is the equiv- 
alent hydrogen column. Both components of the partial coverer 
were multiplied by the cabs model using the same column den- 
sities, (Vh,j and Nh, 2 , respectively. The Gaussian components 
used to model the fluorescent lines were only photo-electrically 
absorbed using a galactic IV^* 1 = 1.7 x 10 22 cm' 2 (Kalberla 
et al., 2005), i.e., they were modeled as they appear outside 
the absorber. 3 This description was chosen because the spectrum 
shows clear Ar and S Ka lines at 2.3 and 2.9 keV, respectively, 
which can not be as strong as observed if they have to pass 
through an absorption column of ~ 10 24 cirT 2 , as is necessary 
to describe the continuum. For a discussion on self-absorption 
and intrinsic strength of the lines, see Sect. 5.2. This spectral 
model, including the same emission lines, can also describe the 
broad-band spectrum of GX 301-2, as seen in recent Suzaku 
data (Suchy et al., 2011). Spectral parameters can be found in 
Table 1 and a plot of the spectrum with the separate components 
and the residuals in Fig. 8. In Table 1 Nh.i gives the absorption 
column which is multiplied by the covering fraction (CF), while 
(Vh ,2 is multiplied by (1 -CF). The normalization and the photon 
index of the power law continuum are given in rows “norm” and 
“photon index F\ respectively. 

The Fe Ka line shows strong residuals when fitted with a sin- 
gle Gaussian. Watanabe et al. (2003) have shown that the Fe Ka 
line shows an extended Compton shoulder (CS) at its lower en- 
ergy side, but the resolution of EPIC-pn does not allow one to 
identify this feature uniquely. Bianchi et al. (2002) have shown 
that the Compton shoulder can be described with a box-shaped 
model. We included such a box in our model, which extends the 
Gaussian shaped Fe Ka line to lower energies and starts at the 
centroid energy of the line. This model gives a very good de- 
scription of the line, as evident in Fig. 8. In Table 1 the row “CS 
Ei” gives the low energy end of the box, while "CS A” gives 
its area. The ratio of the flux between the Fe Ka line and the 
Compton shoulder is in good agreement with the ratio ex- 
pected for photons transmitted through a dense absorption 
column on the order of 10 23 cm" 2 (Matt, 2002). 

2 see also http://puisar.sternwarte.uni-erlangen.de/ 
wilms/research/ tbabs/ 

3 DL value obtained with nh (HEASOFT v6.18). provided by 
HEASARC, at the coordinates of GX 301 -2 


The best fit x~ obtained with this model was 769 with 376 
degrees of freedom (dof) 0/ 2 cd = 2.04). Still snake-like residuals 
between 5.0-5 .5 keV were visible, see panel h) of Fig. 8. The 
Ka line energy of neutral chromium is at 5.41 keV (Kaastra 
& Mewe, 1993; Kortright & Thompson, 2001), and adding a 
Gaussian emission line in this energy range let to an improve- 
ment of the fit to x 2 = 713 with 374 d.o.f. The center of the 
line was found at 5.425"®'* 13 keV in agreement with theory. Still 
residuals remain around 5.6 keV in absorption, see panel c) of 
Fig. 8. These residuals are in the energy range of the plateau be- 
tween the iron lines and their escape peaks in the silicon of the 
detector. The escape peaks and the plateau are clearly visible in 
the model of the Gaussian components of the Fe Ka and Fe K/3 
line in Fig. 8 due to the large flux of these lines. The residuals 
could be very well described with a broad Gaussian absorption 
feature, improving the fit further tox 1 = 488 with 371 d.o.f. and 
called “Abs. Feat.” in Table 1 . Using only this absorption feature, 
and leaving out the Cr Ka line, to model the residuals does not 
adequately describe the spectrum, and gives a^ 2 of only 562 for 
373 d.o.f., see Fig. 8d. 

Finally strong residuals remained around 8.3 keV, very close 
to the energy of the Nickel K/3 line at 8.26 keV (Palmeri et al., 
2008a) . A Gaussian shaped line centered at that energy was used 
to model the residuals and improved the fit clearly. This last ad- 
dition resulted in our best fit model , with x 2 = 448 for 369 d.o.f. 
{X 2 rci = 1.214), of which the residuals are shown in Fig. 8 f. 

To our knowledge, this is the first time that clear evidence for 
a Cr Ka as well as a Ni K/3 emission line was found in the soft X- 
ray spectrum of a HMXB. Cr has been observed in absorption in 
other sou rces before , e.g ., in G RO J 1 655-40 (Miller et al . , 2008) . 
These authors found, however, no evidence for fluorescence lines 
in the spectrum of GRO J 1655-40 and Cr as well as all other Fe- 
like elements were found to be highly ionized. 

Even though the spectrum was only used above 2 keV, as 
all bins below are clearly background dominated, the continuum 
in the energy range between 2-3 keV can not be accurately de- 
scribed by only the partial covering model. A constant DC level 
is needed additionally, to lift the continuum in the lowest energy 
bins. Without this DC component the best fit provides only a^- 2 - 
value of 737 with 370 d.o.f. Values of the DC level are given 
in Table 1 in the row labeled “DC level”. The DC level might 
originate from a low energy response of the large Fe Ka line be- 
low the escape peak, which is not fully taken into account in the 
RMF. A similar effect is known to be present in XIS data from 
Suzaku (Matsumoto et al., 2006). 

4.2. Time resolved spectra 

We have seen in the variability of the light curves, pulsed frac- 
tions, and hardness ratios that the overall spectrum will only 
provide a very inadequate description of the physics of the neu- 
tron star, as many different accretion regimes are mixed together. 
Sacrificing S/N we therefore decided to split the spectrum ac- 
cording to the four regimes I-1V defined in Sect. 3 and indicated 
in Fig. 2. These states show the strongest differences in the light 
curve and hardness ratio and are therefore expected to show dif- 
ferent spectra. 

All spectra were rebinned to a S/N of 12 up to 3 .5 keV, of 36 
up to 6keV, of 60 up to 8.5 keV and of 48 above. Additionally, 
at least 4 channels were combined to account for the real energy 
resolution of the instrument. This binning resulted in different 
amounts of degrees of freedom for each spectrum. To each of 
these spectra we fitted the same model as to the overall spectrum 
and obtained very good values ~ 1 for all four spectra. 
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Fig. 8. a): Phase and time averaged spectrum in the 2-10 keV energy 
range of EPIC-pn. The components of the spectral model are shown 
in different colors, see text for a detailed description, b) Residuals 
with neither an absorption feature nor an Cr K a line around 5.4 keV. 
c) Residuals with a Cr Ka- line, but without an absorption feature, d) 
Residuals with an absorption feature, but without a Cr Ka line, e) 
Residuals without a Ni K/J line./) Best fit residuals. 


The best fit parameters are listed in Table 1 . It can be seen 
that the absorption is highest during part II in which the coun- 
trate drops very low and the pulsations vanish. During that time 
the absorption is increased by ~20% compared to part I. The un- 
derlying power law continuum softens at the same time, result- 
ing in comparable hardness ratios. The unabsorbed flux is also 
highest during part II (see row y, u _“^ eV in Table 1), so that the 
observed drop in count rate can be attributed solely to the very 
high absorption column. 

The width, cr. of the Fe Ka was in all parts around 40 eV, 
larger than the energy resolution of EPIC-pn. If this width 
is due to Doppler shifts from turbulence in the fluorescent 
medium, speeds in excess of lOOOkms -1 would be necessary. 
In GX 301-2. observational evidence and models suggest that 
a wind speed of only 300kms~ ! is reached at the position of 


the neutron star (Kaper et al., 2006; Leahy & Kostka, 2008), 
which can account for significantly lesser broadening of only 
~10eV. During accretion the wind will become turbulent, how- 
ever, and thereby possibly locally increasing the wind speed to 
values high enough to explain the broadening (Mauche et al., 
2008). Turbulent flows in the accretion stream are, however, the- 
oretically not yet well understood. A more likely source of ap- 
parent line broadening comes from the fact that the Fe Ko- is 
not a single line, but a superposition of different K a lines of dif- 
ferently strongly ionized iron (Kallman et al., 2004). This effect 
can account for a width of 20-30 eV of the observed line, as the 
single lines cannot be resolved with the spectral resolution of 
EPIC-pn. Together with broadening in the wind, the observed 
width of Fe K a seems realistic. The Fe K/J line is even broader 
with values around 65 eV. As the Fe K/J complex is not as broad 
as the Fe Ka complex (Kallman et al., 2004), this broadening 
might be mainly due to confusion with the strong Fe K-edge at 
7.1 keV, very close to the line. EPIC-pn’s resolution does not al- 
low to completely separate these two features, but attempts to 
fit the spectra with a narrow Fe K/J line lead to unacceptable 
T'd values - 

The centroid energy of the Fe K a line is changing with time 
and can be as high as 6.45 keV in part III, where the brightest 
flare occurred. This energy corresponds to Fexix or higher ion- 
ized iron, which does not emit a Fe K/J line, as all M-shell elec- 
trons are already stripped from the ion. Still the spectrum of part 
III shows clear evidence for a Fe K/J line at 7.127keV, corre- 
sponding roughly to Fe xm to xiv (Kallman et al., 2004). For the 
Ni Ka and Ni K/J lines a similar discrepancy emerges regard- 
ing the best fit values: the best fit Ni K/3 centroid energy is very 
close to the one of neutral nickel (Palmeri et al., 2008a), while 
the Ni Ko- centroid energy requires at least Ni X or higher. The 
uncertainties of the Ni K/J energy, however, allow also for ion- 
ization states of this order. 

Regarding the flux ratios between the Fe Ka- and Fe K/J we 
obtain strongly varying values between Fe K/J/Fe Ka = 0.186 - 
0.283. The Ni K/J/Ni Ka ratio seems to be less variable with 
values between 0.18 - 0.23, and the individual uncertainties 
are much larger so that they are all consistent with each other 
(see Table 1, rows “Ratio Fe“ and “Ratio Ni“). These values 
are clearly above the theoretical ones for neutral elements from 
Kaastra & Mewe (1993), which give Fe K/J/Fe Ka = 0.125 and 
Ni K/J/Ni K/J = 0.105. Measurements in solid state metals give 
higher numbers of Fe K/J/Fe Ka = 0.132 and Ni K/J/Ni Ka 
= 0.133 (Han & Demir, 2009). The Ni K/J/Ni Ka ratio is 
marginally consistent with all these values within its uncertain- 
ties, but the measured iron ratio is clearly higher. Other measure- 
ments and calculations for neutral iron give figures of the same 
order with Fe K/J/Fe Ka « 0.13 (see, e.g., Palmeri et al., 2003, 
and references therein). 

It is conceivable that the flourescence lines originate 
from the dense absorbing material. If this is the case, self- 
absorption is not negligible (Inoue, 1985 ). For an column den- 
sity of A h a 140 x 10 22 cm -2 the optical depth at the energy 
of the Fe Ka line is only ~ 60% of the depth at the energy of 
the Fe K/J line. If the whole absorption column would act on 
the lines, the intrinsic ratios are drastically reduced to val- 
ues between Fe K/J/Fe Ka = 0.114 - 0 . 173 , much closer to 
the theoretical values. Such a strong effect is, however, un- 
likely, as the fluorescence lines at softer energies from ele- 
ments like S and Ar would vanish completely. The line ra- 
tio also depends strongly on the viewing angle and the ge- 
ometry of the absorber (Molendi et al., 2003 , and references 
therein). Simulations have show n that for certain viewing an- 
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Table 1. Fit parameters for the best fit model for the overall spectrum and the spectra of the four states i-IV, respectively. For details on the model 
see text. Uncertainties are given at the 90% limit. The two lowest values are obtained when using a gainshift kernel, as discussed in the text. 


Parameter 

all 

I 

11 

III 

IV 

Cov. Frac. 

0 935 +0 005 

0 943 +uon 
-0.012 

0 9796" 0000 

0 89~°-° 2 

^■° -0.03 

0.968";“!“ 

An i [ 10 22 atoms crrr 2 l 

141*1 

147 ±6 

175" 

114 + 5 

130"- 2 

Nil 2 [ 10 22 atoms cm -2 ] 

53.7 ± 1.6 

59 ±4 

47 ±2 

53" 

42 :i 

Norm [keV' 1 s' 1 cm' 2 ] 

t | r+0.00 
1,1 -0.06 

i iq+0.18 
1-1 -0.15 

2 2* 03 
^- 0.1 

1 A7+0.09 
1 U ^-0.10 

l.l" 4 

Photo Index F 

0.90 ± 0.02 

0.93 + 0.04 

1.10 + 0.04 

0.78 ± 0.03 

1-055S 

Fe Ka A fph s' 1 cm' 2 ] 

(l.04!»;g)xl0" 2 

(0.82:«")xl0~ 2 

(i.io:“ 3 ) x io~ 2 

(l.27!“;*)xl0' 2 

(0-78!“;®) x 10" 

FeKaE [keV] 
FeKacr[keVl 
CSE, [keV| 

6A22T 4m 

0 A426"” 018 
u.o+zo_ ooois 

6.205 ± 0.009 

6.413 ± 0.003 
0.040 + 0.004 
6.199!“" 

6.407 ± 0.002 
0.040 ± 0.003 

< ioi +0.013 
-0.015 

6.453 ± 0.003 

0 032 +0003 
6-249!^ 

6.4inr 

rXSsr 3 

u -^ J -0.04 

CS A [phs' 1 cm' 2 ] 

(2.0 ± 0.3) x 10' 3 

(2 4”; 4 ) x 10 1 

(2.2!“ 4 ) x 10" 3 

(3.3!“)x 10" 

(l.3!“; 2 )xl0' 3 

Fe K/3A [phs' 1 cm' 2 ] 

(2.32 ± 0.05) x 10' 3 

(1.76 ± 0.08) x 1 O' 3 

(2.04 + 0.07) x 10" 

(3-68!«;J 3 ) x 10-3 

(l .53!g;!3) x 10" 

Fe K/3E [keV] 

7.095 ± 0.002 

7.087 ± 0.005 

7 O70 +0003 

7 12669 + °- 00009 

7.063 ± 0.009 

FeKfScr [keV] 

0.0675"®} 

0.069 ± 0.006 

0.063 ± 0.005 

0.052 + 0.005 

o.o58:“!!! 2 

Ni Ka A [ph s' 1 cm' 2 ] 

10 ~ 3 

(4.4^) x 10' 4 

(0.66";“)xl0" 

(0.69 ± 0.07) x 10" 

(o.42:“ !- 3 ) X 10" 

NiKaElkeV] 
Ni Kacr [keV| 

7.491 ±0.005 
0.036 + ““!? 

7 4S4+ 0 - 005 

/ - 4+ O^_0.011 

<0.04 

7 4772 + 0.004 
-Q.Q0Q 

0 053 +9 - 91 ‘ 
U,UJJ - 0 . 0 1 2 

7 51 5' 0 ' 015 
< 0.03 

7.50 ± 0.03 

0 06 + ° 05 
0.06 

CaKttA [phs -i cm -2 ] 

(0.73 ± 0.06) x 10 " 4 

(0.47 ± 0.09) x 1 0 “ 4 

(0.82 ± 0.09) x 10 " 4 

(1.03 ± 0.18) x 10 " 4 

( 0.8 ± 0 . 2 ) x 10' 4 

CaKaE [keV] 

o-fi + 0.02 
J * / 1 - 0.01 

771 +0.03 
Jm ' *- 0.00 

3 716 +0 004 
0/1 °- 0.01 1 

2 7 Q <+0.015 
j. / JJ_ 0 009 

3 69 +0 02 

-0.03 

ArKaA |phs ' 1 cm' 2 ] 

(3.2 ± 0.3) x 10" 

(1.9 ± 0.5) x 10" 

(3.1 ±0.5) x 10" 

(0.47 ± 0.09) x 10 " 4 

f0.48 + “ !“'l x 10 ' 4 

ArKaE [keV] 

2-973:“;“; 

2.969!“|“ 

2 <HS +0 - 015 
- 0.000 

3.000 ±0.015 

2.96 ± 0.02 

S Ka A [phs " 1 cm' 2 ] 

(3.3 ± 0.4) x 10 ' 5 

(2.1 ± 0.5) x 10 ' 5 

(3.0 ± 0.5) x 10" 

(0.61 ± 0 . 10 ) x 10' 4 

( 2 . 9 r; |) x 10 " 

S Ka E [keV | 

2-336!“;“g 

2.35 ± 0.02 

2.315 + 0.011 

2.348 + 0.012 

2.36 ± 0.06 

Abs. Feat. A [phs ' 1 cm" 2 ] 

(-0.52:“;?“) xlO ' 3 

(-3.3!! 4 )x 10 ' 4 

(-2.3!“) x 10 ' 4 

(-1.5 + 0.5) x 10" 

(- 0 . 8 " 3 ) x 10 " 

Abs. Feat.E [keV] 
Abs. Feat.cr [keV] 

5.61 ±0.03 
0.26 ± 0.03 

c <;7+0.06 

o-is+l 

5 67 +0 03 
' - 0.02 

0 12 + ^ 11 
y- lz -o.o5 

5.7°!“! 

0-45!“;® 

5 56 +ai9 

~’-~’ u -Q.08 

0 32 +9 - 16 

U,JA -0.I2 

CrKaA [phs -1 cm' 2 ] 

(1.5 ± 0.2) x 10 ' 4 

(1.2 ± 0.5) x 10 ' 4 

(o.9!“- 2 ) x 10 ' 4 

(l.9!“)x 10 ' 4 

(1.9 + 0.7) x 10 ' 4 

Cr Ka E [keV] 

5 430 + °° 15 
J,HJU - 0 . 0 00 

c 45 + 0 .M 

5.42 ± 0.03 

5.44 ± 0.03 

3 A A +0-02 

NiKMlphs-'cm' 2 ] 

(1.2 ± 0.3) x 10 ' 4 

(0.8 ± 0.5) x 10 ' 4 

(1.2 + 0.4) x 10 ' 4 

(1.9 ± 0.8) x 10 " 4 

( 0 . 8 !“; 9 )xl 0- 4 

NiK/3E [keV] 

8-26"“ 

0 00+0.06 
®-^_0.08 

8 - 22 !“S 

8.27!“;“ 

8 . 2 "j 

DC Level [keV ' 1 s ' 1 cm' 2 ] 

(0.48 ± 0.04) x 10 ' 4 

(3.9 ± 0.6) x 10 ' 5 

(3.4 ± 0.6) x 10" 

(0.74 ± 0.11) x 10 ' 4 

(0.4" 3 ) x 10 ' 4 

Ratio Fe 

0 223 +a007 
V-^_0.006 

0 214+ 0014 

-0.013 

0 -i 86 ;“™ 

A 700 + 0.020 
0.283_ 0016 

0 . 20 !“" 

Ratio Ni 

0.20 ± 0.05 

0.18 ± 0.11 

0.18 + 0.06 

0.23 + 0.14 

0.2 ± 0.2 

Fe Ions 

- 

XII-XIII 

VI-X1 

>X1V 

IX-X1 

^ 2 /d.o.f(AL ) 
? 2 - lOkev [erg sec ' 1 cm"-] 

448.02/369(1.214) 

1 13.69/122 (0.932) 

136.13/120 (1.134) 

129.76/131 (0.991) 

91.72/90(1.019) 

1.043 x 10 ' 9 

8.863 x 10 ' 10 

7.127 x 10 ' 10 

1.852 x 10" 

8.048 x 10 ' 10 

TTfmLv Lerg sec ' 1 cm' 2 ] 

1.617 x 10" 

1.624 x 10 ' 8 

2.151 x 10 ' 8 

1.797 x 10" 

1.350 x 10" 

Gainshift 

1.00396 ±0.00017 

1 >(XXU 

i.uwiu _ 0(XX)2 

1.0082 + 0.0003 

1 00 19 + 0.0006 

1,UUA -0.0005 

X 2 I d-o.f (xD 

468.087375(1.248) 

129.48/128 (1.012) 

172.64/126 (1370) 

133.93/137 (0.978) 

101.12/96(1.053) 


gles around 45° of a dense slab absorber in transmission the 
measured ratios are reduced, while softer fluorescence lines, 
especially the Cr K a line, are still clearly detected (Barragan 
et al., 2011, in prep.). A detailed investigation will be pre- 
sented in a forthcoming publication. For further discussion 
of the absorption of the fluorescence lines see also Sect. 5.2. 

The line energies as well as the flux ratio is expected to 
change with the ionization state of iron (Palmeri et al., 2003). 
Allowing for a mixture of different ionizations during each part 
and using the energies and yields from Mendoza et al . (2004), we 
found that ions between Fe vi to Fe xiv can explain the observed 
energies and ratios (see Table 1 , row ”Fe lons“). Here we ne- 
glected self-absorption, as we can not confine its effect on the 
ratios. Additionally the Fe K/J line is superimposed on the iron 
K-edge and the Gaussian model component can be artificially 
increased due to a correlation with the edge. The real ioniza- 
tion states are therefore likely to be lower, but changes are 
still clearly present. We find an indication that the ionization 
state is correlated with the luminosity of the source but that 
overall the ionization state stays around intermediate values. 


It is furthermore possible that the energy calibration of XMM- 
Newton is slightly off and that a small shift in energy depending 
on the count-rate is present, leading to an apparent shift of the 
observed line energies. The most likely effect causing such a be- 
havior is either a change in the charge transfer inefficiency (CTI) 
of the CCD or a gainshift due to X-ray loading. The CTI is de- 
creasing with higher photon flux when impurities in the CCD are 
filled up with electrons (Striider et al., 1997), while X-ray load- 
ing causes more electrons than expected. If these changes are not 
taken into account correctly, they lead to a shift of measured en- 
ergies depending on the source flux. The expected uncertainties 
of the calibration of these effects in EPIC-pn are below ±0.5% 
(Guainazzi, 201 1 , priv. comm.). 

To investigate the possibility of such a shift we froze the 
line energies of all measured elements to the values of their 
respective neutral line energy, with values taken from Palmeri 
et al. (2008b) for S, Ar, and Ca, Kaastra & Mewe (1993) for 
Cr, Kallman et al. (2004) for Fe, and Palmeri et al. (2008a) for 
Ni. We then applied a gainshift kernel, shifting the model counts 
spectrum energy grid according to E ncvt (c) = £ 0 ]d(c)/gainshift. 
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0 10000 20000 30000 40000 

Time [sec] 

Fig. 9. Evolution of the spectral parameters of individual pulses, a) 
Light curve between 0.5-10 keV , b) Fe Ko- flux (including the Compton 
shoulder), c) Fe K/3 flux, d) Ni Ko flux, e) ratio between the Compton 
shoulder and the Gaussian shaped Fe K a line, and f) ratio between the 
Fe Ko and the Fe K/3 flux. The fluxes of the lines are given in pho- 
tons s' 1 cm -2 x 10~ 3 . 


The obtained best fits were only marginally worse than with free 
energies of the lines and the shifts were on the order of 1 .001- 
1 .008 (0.1-0 .8%), see bottom lines of Table 1 . These values are 
only marginally higher than expected from calibration uncertain- 
ties but also include unknown source intrinsic energy shifts. The 
continuum parameters of the spectral model did not change sig- 
nificantly due to the gainshfit, but the fluxes of the fluorescence 
lines were slightly altered to values marginally consistent with 
the ones shown in Table 1 . Especially the Fe K/J line flux was de- 
creased, possibly due to a better seperation from the iron K-edge. 
This effect influenced also the Fe K/3/Fe Ka ratio by decreasing 
it slightly to values <0.2 at all times, which is closer to the values 
expected for neutral iron. No significant change with luminosity 
could be observed. 

Using the gainshift kernel we can explain the measured line 
energy shifts stemming purely from detector effects while the 
physical condition of the plasma in GX 301-2 does not change 
significantly during our observation. 

4.3. Pulse to pulse spectroscopy 

In an approach to understand the time variability of the spec- 
tral parameters, especially the iron line complex, in more detail. 


0 10000 20000 30000 40000 

Time [sec] 

Fig. 10. Time series of the spectral parameters of the spectra of indi- 
vidual pulses, a) Light curve between 0.5-10 keV, b) continuum norm 
in photons keV" 1 s' 1 cm -2 , c) main absorption column in 10 22 cm~ 2 , d) 
covering fraction, e) gainshift slope, and t * le ' 5est 


we extracted spectra for individual pulses, leading to 69 spectra 
with 685.0sec exposure time each. The resulting model param- 
eters are presented in Figs. 9 and 10. The data are grouped by 
color into 9 parts labeled A-I to facilitate orientation. We ap- 
plied a reduced version of the time averaged spectral model and 
included only the Fe Ka line with its Compton shoulder, the 
Fe K/3 , and the Ni Kcr line. All other fluorescent lines could not 
be identified in these short exposure spectra. Additionally we 
had to fix the photon index T = 0.90 and the secondary absorp- 
tion column jVh ,2 = 5.37 x 10 23 cm -2 because neither could be 
constrained well with the data. When leaving the centroid energy 
of the emission lines free, we found a strong correlation between 
these energies and the average count rate, as already indicated in 
the previous Section. This correlation is clearly stronger than ex- 
pected from ionization (Kallman et al., 2004) and is therefore 
likely a detector effect. As described in the previous Section, 
these calibration uncertainties can be modeled with a gainshift 
kernel. We also checked the spectra extracted using epfast, but 
they showed the same systematic effects. 

With fixed energies for the fluorescent lines and for the iron 
K-edge the use of a gainshift kernel leads to a very good de- 
scription of the individual spectra. As expected the gainshift is 
directly correlated to the incident countrate, see Fig. 10c. 

The flux of the Fe K a line and its Compton shoulder is 
mostly, but not always, correlated with the count rate in the 0.5- 
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lOkeV energy band, as seen in Fig. 9. Especially in the second, 
main flare (part G) the iron line flux starts to rise with the light 
curve but breaks down around the maximum of the flare and 
stays low for the remainder of the observation. A better corre- 
lation is observed with the normalization of the power law con- 
tinuum (Fig. 10). The Fe K/3 and the Ni Ka lines show a very 
similar behavior and are strongly correlated with the Fe Ka line 
flux, with a Pearson’s correlation coefficient of 0.92 and 0.89, re- 
spectively. As discussed in the previous Section, the Fe Ka and 
Fe K/3 lines required a finite width <r, with an average value of 
40 eV and 71 eV, respectively. The widths did not show a signif- 
icant correlation with time or flux. The width of the Ni Ka was 
frozen to 1 0“ 4 keV. 

The absorption column is also strongly changing over the 
course of the observation, varying by more than 50% (see 
Fig. 10c). It is increasing during the first flare in part C and stays 
high until the onset of the second, main flare (parts F and G), 
where a strong drop is registered. After the flare it is increas- 
ing again gradually. The covering fraction is very high through- 
out the observation, with an average of ~92%, reducing the sec- 
ondary absorption column to almost negligible values. It is only 
during the bright flare in part G where the covering fraction sud- 
denly drops as low as 74%, around the time when the Fe Ka line 
flux dips. 

A typical spectrum with one pulse exposure time is shown in 
Fig. 11. The Fe Ka and Fe K/3 lines as well as the iron K-edge, 
are clearly seen. The Ni Ka line is also visible on close inspec- 
tion, but below 3 keV no significant excess over the background 
can be measured anymore due to the huge absorption column. 
The Ca, Ar, S, and Cr Ka lines and the Ni K/3 line are not signif- 
icantly detected. 

Assuming that the absorbing matter is the same matter re- 
sponsible for the fluorescent iron emission and that it is dis- 
tributed spherically symmetric around the X-ray source, the 
equivalent width of the Fe Ka line is proportional to the over- 
all absorbing column (Inoue, 1985). Furthermore, self-shielding 
should be negligible and the matter is assumed to be cold and 
neutral. We can see this behavior during most of the XMM- 
iVewton-observation, even at very large Ah values, see Fig. 12. 
The average Pearson’s correlation coefficient over all parts ex- 
cept part D is 0.91 . During part D, corresponding to the off-state 
seen in the light curve (Figs. 2 and 5), this correlation is not vis- 
ible. Instead indication for an anti-correlation with a correlation 
coefficient of -0.67 exists. During that state the equivalent width 
of the Fe Ka line reaches its highest values, whereas the Ah 
stays in a regime around 1 .8xl0 24 cm -2 , lower than the maximal 
measured value. The overall Ah value is quite high, as expected 
from the analysis of region II from Sect. 4.2 (corresponding to 
parts D, E, and F) which showed the highest overall absorption 
value. But the off-state is happening during part D only, a part 
which does not stand out for its increased absorption column 
when compared to parts E and F. This invariance of the Ah value 
clearly shows that the off-state is not due to a increased absorp- 
tion, which was already indicated by the softening during that 
state (see Fig. 2). 

5. Discussion & Conclusion 

We have presented a detailed timing and spectral analysis of the 
HMXB GX 301 -2 during its pre-periastron flare using the EPIC- 
pn camera of XMM-Newton . The 2-1 0 ke V spectrum can be de- 
scribed with a partially covered, strongly absorbed power law 
model enhanced by several Gaussian emission lines from S, Ar, 
Ca, Fe, and Ni. We found evidence that the Ka line from Cr is 



Fig. 11. Typical spectrum with 685.0 sec exposure time, i.e., one pulse. 
The upper panel shows the data and the best fit spectrum, the lower 
panel shows the respective residuals. For details about the spectral 
model see text. 



JV H ,i [10 22 H-atomscm~ 2 ] 


Fig. 12. Equivalent width of the Fe Ka line as function of the main 
absorption column Ah,i ■ Colors are the same as in Fig. 10 and are indi- 
cating the time. 


also present as well as the K/3 line from Ni. To our knowledge 
this is the first detection of these lines in a HMXB. In order to ac- 
count for the different states of activity visible in the light curve 
we separated the data according to these states, labeling them 
I— IV. During the sub-flares (parts I and III) regular pulsations 
with P = 685.0sec were clearly detected and a double peaked 
pulse profile is visible. During the off-state II these pulsations 
vanished almost completely and the countrate dropped to very 
low values of ~10% of the peak count rate. XMM-Newton still 
detected the source clearly. By separating the data into individ- 
ual spectra for every rotation of the neutron star, we were able 
to perform a detailed spectral analysis of the off-state. We found 
that although the absorbing column is very large, it is not dras- 
tically increased during that state, compared to values measured 
shortly before and after it. Instead a softening of the spectra is 
observed in which the Fe Ka line is even more important than at 
other times. 
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We have shown that the absorption column is larger than 
10 24 cm" 2 at all times during the observation. GX 301-2 is 
not the only source showing absorption columns of this magni- 
tude. INTEGRAL Jor example, discovered more than 10 highly 
absorbed HMXB, which might form a new class of objects 
(Kuulkers, 2005; Walter et al., 2006). Among them are sources 
like IGRJ 163 18-4848 (Matt & Guainazzi,2003; Barragan et al., 
2009) or IG R J 1 7252-36 1 6 (Manousakis & Walter, 2011) which 
show very prominent Fe Ka and Fe K/3 lines, as well as Ni K a 
lines. The flux of these sources is, however, usually lower than 
that of GX 301-2, decreasing the S/N drastically below 4keV 
and thereby not allowing for studies of softer fluorescence lines. 
Results from the iron line energies and Fe Ka/K/3 ratios also in- 
dicate only mildly ionized iron. GX 301-2 could therefore be 
a brighter example of these sources and help to understand this 
class of HMXB. 

5. 1 . Fluorescence region size 

To estimate the distance between the X-ray producing region and 
the fluorescent medium, we extracted light curves with 1 sec res- 
olution for the iron line energy band between 6.3-6.5keV and 
the harder continuum band between 7.3-8.5keV. The countrate 
between 6.3-6 .5 keV is completely dominated by the Fe K a line 
and any variability of it must be due to the reaction of the flu- 
orescent line to changes in the continuum flux. We compared 
both light curves using Fourier transformation as well as cross- 
correlation techniques and could not find a significant time de- 
lay between both curves in the range of 2-5000 sec. This is in 
agreement with the size of the fluorescence region being smaller 
than 10 10 cm, i.e.,0.3 lt-sec, as given by Endo et al. (2002). The 
count rate in the narrow energy bands was not high enough to 
achieve a time resolution sufficient to see a delay of the fluores- 
cence lines on that order. For IGR J 17252 — 3616 Manousakis & 
Walter (2011) also found that the emission line region is com- 
pact and not larger than the accretion radius. However, during 
the off-state in part D in Fig. 9 the iron line is still clearly visible, 
although the luminosity of the continuum has significantly de- 
creased. This leads to the largest equivalent widths for the Fe K a 
line measured in the whole observation (Fig. 12) and seems to 
indicate that the fluorescence line is partly reacting rather slowly 
on changes in the continuum. A possible explanation is that the 
weaker part of the partial coverer is located further out, and 
thereby reacting much slower to any changes in the X-ray contin- 
uum flux. The X-ray flux was low for ~4500sec and we can as- 
sume this scale as an upper limit of the light travel time to the ab- 
sorber. On a direct path, this would place the absorbing medium 
almost 10 times the average orbital diameter of the system away. 
Through scattering in the dense medium, with an optical depth 
r = Nh.i x ctt = 1, the effective light-path is increased, so that 
a size of 2000 lt-sec or less is plausible. A similar geometry was 
already proposed by Mukherjee & Paul (2004). Through scatter- 
ing in the inner absorber the outer absorber receives a much less 
pulsating X-ray flux and transforms it into fluorescence lines, 
which are consequently also pulsating much less. This effect can 
explain the reduced pulsed fraction in the Fe K a and Fe K fi line 
band (Fig. 7). A secondary absorber further out can also explain 
the dip in N H and covering fraction during the large flare in part 
G (Fig. 10). Here the inner absorber is almost completely ion- 
ized and transparent to soft X-rays. The visible absorption is 
then primarily due to the still cold outer absorber. The present 
data do not allow one to distinguish between the lines formed in 
the different absorbers, and as they are both variable it is hard to 
disentangle the effects. The indications we obtain from the data 


do, however, support this absorption geometry. If both columns 
of the partial covering model would be located close to the neu- 
tron star, where the strong column describes the clumpy stellar 
wind and the other the smooth part of the wind, neither the large 
equivalent width nor the strong dip in covering fraction during 
the large flare is easily explainable. The measured widths of the 
Fe Ka and Fe K/3 lines as well as their ratio (Sec. 4.2) can also 
be understood in this geometry as a superposition of cold and 
neutral iron in the outer absorber and higher ionized iron in the 
inner one. As stated before, the spectral resolution of F.PlC-pn 
does not allow to separate different iron lines, but we would ex- 
pect to see an only slowly variable Fe Ka line of neutral iron and 
many variable lines from the different ionized species from the 
inner absorber. 


5.2. Fluorescence and abundances 

The 2-1 0 keV spectrum of GX 30 1 -2 shows evidence for Fe Ka, 
Fe K/3 , Ni Ka, Ni K/3, S Ka, CaKa, ArKa, and CrKa fluores- 
cent emission lines. To check consistency with older data, es- 
pecially regarding the CrKa line, we applied our model to the 
Chandra HETG data presented by (Watanabe et al., 2003) and 
taken during the pre-periastron flare in January 2001 . The model 
describes the data very well in its overall shape. Neither a Cr Ka 
nor a Ni K/3 line is significantly detected, but lines with the same 
flux and width as in the XMM-Newton spectrum are also not 
ruled out by the data and could vanish due to the worse S/N of 
the Chandra data. 

To describe the XMM-Newton spectrum we used a model in 
which the fluorescent lines were only absorbed by the galactic 
absorption column of N = 1.7 x 10" cm~ 2 to accommodate 
for the observed fluxes of the Ka lines of S and Ar, see Sect. 4.1 . 
The flux of the lines indicates, however, that they originate in 
a dense absorber and that thus self-absorption must be taken 
into account. Additionally timing analysis has shown that the 
main producing region of the iron line must be very close to the 
X-ray source and it is likely that this is true for the other lines 
as well. It is thus conceivable that the photons emitted from the 
fluorescent ions are also absorbed by material relatively close to 
the neutron star with a column density in excess of A* al , e.g., by 
the stellar wind of the companion. 

In order to estimate how strongly these lines are absorbed 
we performed the following simple sanity check, assuming so- 
lar abundances for all elements (Wilms et al., 2000): for each 
element, where a fluorescent line was observed, we calculated 
how many photons were absorbed by this element alone above 
its respective K-shell , applying the cross-sections by Verner et al . 
(1996). Using the fluorescence yields for near neutral ions from 
Kaastra & Mewe (1993) we calculated how many of these pho- 
tons were reemitted in the respective Ka and K/3 lines. As we 
do not know the geometry of the fluorescent medium and hence 
can not tell how many of these reemitted photons reach us, we 
calculated the fluxes of all lines relative to Fe Ka. We compared 
these relative fluxes to the fluxes as measured by the Gaussian 
lines in our model. The relative fluxes of the Gaussian lines are 
strongly correlated with the absorption column applied, espe- 
cially for the soft S and Ar lines. By changing the absorption 
column we can find the best agreement between measurement 
and calculation. We find that this is achieved with an overall ab- 
sorption column around 20 x 10 22 cm -2 , i.e., a factor of ~6 less 
than the continuum absorption and consistent with the above as- 
sumptions. Systematic uncertainties outweigh any statistical un- 
certainties in this analysis and are estimated to be on the order of 
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±5 x 10 22 cm“ 2 . A drastic improvement on this estimate can 
be achieved using Monte Carlo simulations of different ab- 
sorber geometries, but this investigation is beyond the scope 
of this paper (see Barragan et al., 2011, in prep.). 

5.3. The chromium line 

We have reported the detection of an emission line at 5.43 keV, 
compatible with the Ka line energy of neutral chromium. This 
line is weak and only visible because of the large absorption col- 
umn of GX 301-2 making it bright enough to be visible over the 
continuum. The highly absorbed IGR sources described above 
are typically not bright enough to provide enough S/N to allow 
searches for a Cr line. In Supernova Remnants (SNR) evidence 
for fluorescent chromium K lines was found, most notably in 
W49B (Hwang et al ., 2000; Miceli et al ., 2006) but also in Tycho, 
Cas A, and Kepler (Yang et al., 2009). Cr is a decay product 
of Fe in supernovae, and thus a high iron abundance should be 
accompanied by a high chromium abundance. The abundance 
ratio of Cr to Fe is ~\% for solar mass distribution. In super- 
novae the ratio depends on the explosion model and can be 2% 
or larger for delayed detonation models (see,e.g., Iwamoto et al., 
1999). The SNRs investigated by Yang et al. (2009) show differ- 
ent abundance ratios between 0.63% - 3 .6%, leading the authors 
to suspect that they originate from different type of supernovae 
explosions. Using the method described in the previous Section 
to calculate the relative strengths and therewith the relative abun- 
dances of Cr to Fe we obtain a value of 1.4%, very well within 
the expected values. Taking only the ratio of the measured pho- 
tons in the CrKa and Fe Ka line, we obtain higher values of 
~6.4%, but this value neglects the fact that the K-edge for Cr is 
lower than the one of Fe and that consequently more photons are 
able to ionize Cr. With our measurement we can not infer possi- 
ble explosion models or production channels of the GX 301-2 
system. The measured relative line flux does, however, agree 
with predicted fluorescence yields and element abundances. 

5.4. Accretion during the off -state 

The off-state is perhaps the single most puzzling feature of the 
XMM-Newton data. A similar off-state was only recently dis- 
covered by Gogii§ et al. (2011) in RXTE data. RXTE’s spec- 
tral resolution does not allow for a detailed analysis of the iron 
line region but provides a larger broad band coverage compared 
to XMM-Newton. Gogtis et al. (201 1) fitted individual spectra 
with exposure times of one pulse period each, keeping the iron 
line fixed and letting the photon index of the power law con- 
tinuum vary. They found that the photon index is drastically in- 
creased during the dip, making the overall spectrum much softer. 
It thus seems very likely that the accretion was switched off dur- 
ing that time and that only weak Comptonization in the accre- 
tion columns was possible, while most of the X-ray flux origi- 
nated from the black body on the neutron star’s surface. This re- 
sult is in accordance with a theory developed to explain similar 
off-state in Vela X-l (Kreykenbohm et al., 2008, and references 
therein). According to this theory, accretion ceases during the 
off-state and the source enters the propeller-regime, i.e., a regime 
in which the Alvfen radius is larger than the co-rotating radius 
(Lamb et al., 1973). This transition is triggered by a strong de- 
crease in the mass accretion rate, possibly due to a strongly struc- 
tured or clumped wind (Oskinova et al., 2007; Kreykenbohm 
et al., 2008; Fiirst et al., 2010). Even though GX 301-2 is not 
directly accreting from the wind in the pre-periastron flare but 


rather from an accretion stream (Leahy & Kostka, 2008), strong 
density variations are evident through the flaring behavior of the 
light-curve (Fig. 2). Other observational evidence for a clumpy, 
wind around Wray 977 was put forward by Mukherjee & Paul 
(2004) who studied in detail the time variability of the absorp- 
tion column on short and long timescales. 

All our results indicate that the observed off-state was caused 
by a cessation of accretion and not by an optically-thick clump 
moving through the line of sight. A clump in the line of sight 
would also lead to drastically reduced pulsations and flux, as 
the X-ray flux is scattered in the absorber and thereby smearing 
out the pulsations (Makino et al., 1985). In the spectral analy- 
sis, however, we would expect that the spectrum gets harder, as 
softer X-rays are affected stronger by absorption than harder X- 
rays. We would also expect to see a drastic increase of the ab- 
sorption column. Both are not the case, as seen in Fig. 10. 

We investigated the overall behavior of the pulsed fraction as 
function of Nh, and found almost no correlation (Fig. 13), con- 
tradictory to the results by Makino et al. (1985). Scattering in the 
absorption medium seems to be influencing the pulsed flux only 
marginally, so that the complete disappearance of pulses in part 
D is very unlikely to be caused by larger /Vh values. The variation 
in the pulsed fraction seems to originate from different accretion 
rates, as expected from a clumpy accretion medium. This influ- 
ence is also supported by the variability of the individual pulse 
profiles, as shown in Fig. 3. We do not see a weaker pulsed flux 
with essentially the same shape causing the lower pulsed frac- 
tion, but rather a clear change in the overall pulse profile caused 
by a different accretion rate. 

From results at harder X-rays (Evangelista et al., 2010) a de- 
crease of the pulsed fraction with decreased luminosity would 
be expected. In a simple toy model this can be explained by a 
change in altitude of the X-ray producing region in the accretion 
column (Lutovinov & Tsygankov, 2009). Such a change would 
also be visible in a change of the energy of the CRSF (Suchy 
et al., 2011; Doroshenko et al., 2010). We did not find a corre- 
lation between the pulsed flux and the luminosity in the XMM- 
Newton data. The pulsed fraction during the pre-periastron flare 
is the lowest over the whole orbit and is also much lower in the 
soft than in the hard X-rays (Evangelista et al., 2010). It is thus 
possible that a change in the accretion column might not be visi- 
ble in the pulsed fraction of the 2-10keV band. 

Assuming that the off-state was caused by a cessation of ac- 
cretion, possibly because the neutron star entered the propeller 
regime, the remaining X-ray flux is mainly coming from the 
still fluorescent absorbing medium around the neutron star, as 
described in Sec. 5.1. 

As shown in Fig. 12 during the normal states the equivalent 
width of the Fe K a is correlated with the Nhj value. This correla- 
tion is explained by the assumption that the fluorescent medium 
is the same as the absorbing medium. As a result, a higher ab- 
sorption column leads also to a brighter fluorescence line as 
more atoms in the absorber are ionized by the X-ray continuum. 
During the off-state, the X-ray continuum is drastically reduced 
leading to a reduced Fe Kir line from the main absorber in the 
line of sight. The off-state indicates a very clumpy medium, so 
that the absorption in the line of sight is not necessarily an in- 
dicator for the overall density of the medium surrounding the 
X-ray source. The observed break-down of the correlation be- 
tween equivalent width of the Fe K a line and the /Vh,i value can 
thus be explained by a complex structure of the medium around 
the neutron star. For example, fluorescence from dense medium 
behind the neutron star could arrive simultaneously with the re- 
duced continuum flux at the observe due to different light-travel 
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Fig. 13. Pulsed fraction as function of jV Hit . Error bars on the pulsed 
fraction are arbitrary. Color code is the same as in Fig. 10. 


times. The remaining continuum flux during the off-state could 
also be due to Compton reflection from a dense medium around 
the neutron star, similar to the scenario proposed by Rea et al. 
(2005) to explain a long low-intensity state in the Low Mass 
X-ray Binary (LMXB) GX 4+1. These authors analysed broad 
band BeppoSAX spectra and found that they could be described 
by a cut-off power law plus a reflection component, the latter 
becoming more important during the low intensity state. In this 
state also the Fe Ka line reached the hightes equivalent widths, 
similar to GX 301-2. Rea et al. (2005) found that although the 
pulsations vanished below 7 keV they could still be seen in the 
harder X-rays, an effect we can not investigate with the XMM- 
Newton data but which is not evident in the RXTE / PCA data 
presented by Gogtij et al. (201 1). If the reflecting surface is not 
as well defined as in the geometry proposed forGX 1+4, all pul- 
sations can be smeared out over a broad energy range, so that 
the reflected continuum flux is also unpulsed when escaping the 
medium around the system. More observations in a broader en- 
ergy range are necessary to understand the geometry further. 

5.5. Outlook 

Off-state are seen in an increasingly number of sources, like 
Vela X-l and GX 1+4 and can provide interesting insight in the 
structure of the accreted material as well as the configuration 
of the accretion and absorption region. For GX 301-2 we pro- 
pose that two absorbers exist, a very dense but dumpy primary 
one close to the neutron star, and a weaker and smoother sec- 
ondary one farther away. Fluorescence lines from the secondary 
absorber are much weaker then the ones from the primary and 
only visible during the off-state. High-resolution spectra of the 
off-states from Chandra HETG or future missions will help to 
investigate them in more detail by allowing for better measure- 
ments of the ionization state and the geometry of the absorber. 
The discovery of a Cr Ka line in the spectrum of GX 301-2 
leads to interesting new questions about the abundances of the 
iron family atoms, including Cr and Ni. With higher effective 
areas in future missions, more sources might be found show- 


ing a Cr Ka line and in comparing these findings to GX 301-2 
more insight into the production of HMXB and the preceding 
supernova can be achieved. With broad-band spectroscopy, as al- 
ready possible with the Suzaku observatory, the underlying con- 
tinuum can be constrained simultaneously in a better way than 
with XMM-Newton oder Chandra only. Our collaboration is an- 
alyzing Suzaku data of GX 301-2 taken during orbital phases 
outside the pre-periastron flare (Suchy et al., 201 1). These obser- 
vations provide insight into the behavior of the accretion column 
through exact measurements of the variability of the CRSF en- 
ergy. All this information combined will complete the picture 
of GX 301-2 and possibly also of the other highly absorbed 
sources. 
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